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We study the effect of resonances associated with complex molecular interaction of Rydberg atoms 
on Rydberg blockade. We show that densely-spaced molecular potentials between doubly-excited 
atomic pairs become unavoidably resonant with the optical excitation at short interatomic separa¬ 
tions. Such molecular resonances limit the coherent control of individual excitations in Rydberg 
blockade. As an illustration, we compute the molecular interaction potentials of Rb atoms near 
the 100s states asymptote to characterize such detrimental molecular resonances, determine the 
resonant loss rate to molecules and inhomogeneous light shifts. Techniques to avoid the undesired 
effect of molecular resonances are discussed. 

PACS numbers: 42.50.Dv, 03.67.-a, 32.80.Rm, 78.67.-n 


Rydberg blockade W has recently emerged as a 
promising method for creating and manipulating quan¬ 
tum states of light and matter in applications ranging 
from quantum information processing to quantum 

nonlinear optics nm. The key idea is that strong inter¬ 
action between Rydberg atoms can be used, under cer¬ 
tain conditions, to block the states with more than one 
excited atoms. Multiple Rydberg excitations are sup¬ 
pressed due to level shifts caused by strong long-range 
interactions between Rydberg atoms. This mechanism 
enables performing quantum logic operations between 
atom pairs and manipulate collective many-body states 
of iV-atom ensemble Q] . Such collective states efficiently 
couple to laser fields with the coupling enhanced by a fac¬ 
tor of \fN, see experiments mm- While a number of 
advanced protocols involving Rydberg blockade is being 
explored, an outstanding challenge is to identify and re¬ 
alize conditions for high-fidelity atomic and optical state 
control via Rydberg blockade. 

Here we investigate the effect of molecular resonance 
on quantum state manipulation via Rydberg blockade. 
We demonstrate that the very same interactions that 
cause the level shifts required for blockade also have 
detrimental effects due to a large state density (num¬ 
ber of levels per energy interval) of Rydberg states re¬ 
sulting in a plethora of closely-spaced molecular poten¬ 
tials. Some of these potentials may become, at specific 
interatomic separation, resonant with the driving field 
causing excitations to unwanted doubly-excited Rydberg 
states. While this mechanism was qualitatively pointed 
out mm , detailed understanding of effects of molecu¬ 
lar resonances on collective state manipulation is impor¬ 
tant for high fidelity quantum states control. This is chal¬ 
lenging partially due to the overwhelming complexity of 
molecular potentials especially at small internuclear sep¬ 
arations m ■ Below we demonstrate that the cumulative 
effect on the Rydberg blockade is caused by molecular 


resonances at large interatomic distances where reliable 
theoretical predictions can be made. We derive and com¬ 
pute the rates of resonant conversions to diatoms, show 
that collective qubit rotations are damped, and compute 
the “leakage” and inhomogenous frequency shifts r due 
to diatom conversion. Finally, we discuss techniques to 
suppress the deleterious molecular resonance effects. 
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FIG. 1: Selected S 9 molecular potentials in the 1 GHz win¬ 
dow centered about the 100s + 100s dissociation limit (placed 
at zero energy) for two Rb atoms. The potentials are marked 
by their double-atom dissociation limits at large internuclear 
separations R. Highlighted blockading interaction is the in¬ 
teraction that tunes a pair of 100s Ry atoms away from the 
resonance with driving laser field. The position R x of the 
outer-most resonant molecular potential crossing is marked 
with a circle. Properties of molecular resonances are com¬ 
piled in Table [I] 

We start by computing molecular potentials for two in¬ 
teracting Rb Rydberg atoms by a direct diagonalization 
of the long-range dipole-dipole molecular Hamiltonian. 
On a large energy scale, we find a “spaghetti” of densely 
packed curves exhibiting intricate avoided crossing pat¬ 
terns. The region that is relevant to our discussion is 
centered around the nominally blockaded Rydberg lev- 
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els. As an illustration, we take |r) = 1100s). Consid¬ 
ering that the typical excitation Rabi frequency flo is 
~ 1 MHz we zoom onto a 1 GHz window (Fig. |T|) cen¬ 
tered about the 100s + 100s dissociation limit. In this 
figure the potential that at large R asymptotes to two 
100s atoms is the blockading van der Waals interaction. 
However, we also find several potential curves that at 
short R cross zero energy corresponding to a resonance 
with the laser field. As a result, atoms can be promoted 
into an undesired molecular state corresponding to two 
Rydberg atoms. Properties of these resonant crossings 
are compiled in Table [TJ The outermost crossing with 
the most substantial laser coupling is at R x ~ 6.2 /im. 
Since this value is larger than the average interatomic 
separation for typical experimental number densities m 
Pd ~ 10 11 — 10 12 cm 3 , one may find a fraction of atomic 
pairs inside the volume enclosing molecular resonance 
region. Note that the potential curves were computed 
in the basis of atomic orbitals with orbital angular mo¬ 
menta up to f? max = 2. Increasing f max and adding atomic 
orbitals to the computational basis breed new resonant 
crossings, as the system becomes increasingly chaotic at 
smaller R due to stronger inter-channel couplings and 
thus larger number of avoided crossings. Even our outer¬ 
most resonance can be superseded by crossings at larger 
R, but with suppressed laser couplings. However, the 
parameters of the outermost crossing in Fig. [T] are stable 
with respect to the basis variation. As we demonstrate 
below namely this outermost crossing predominantly af¬ 
fects the dynamics of collective excitations thereby miti¬ 
gating challenges of reliably computing full-scale Ry-Ry 
interaction potentials. 
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TABLE I: Molecular resonance shell properties for the 
Rb 100s blockaded state (see Fig. [I]). R x and A R x are 
the shell radii and widths, are the fractional molecular 

Rabi frequencies, and 7 m are molecular loss 

rates. A R x and are evaluated for atomic Rabi frequency 
Qo = (27 t) x 0.1 MHz and number density pd. = 10 12 cm -3 . 
Positions and the number of resonances for smaller R are 
computational-basis dependent, however, their contribution 
to the total rate is strongly suppressed. 

Despite the complexity of molecular potentials, the po¬ 
sition R x of the outermost resonance can be estimated 
as follows. Suppose the ns + ns state is our |rr) “block¬ 


aded” state. The nearest-energy n's + n"s state with 
n! Ri n" ~ n is the (n — l)s + (n + l)s state and at large R 
it lies below the resonance by S m ~ —3 n~ 4 . Further, the 
molecular potential correlating to the (n — l)s + (n + l)s 
atoms, behaves as n^c^/R 3 due to the repulsion from 
the p + p state below, where C 3 ~ 1. Thereby, U(R) ~ 
—3n -4 + n^ds/R 3 and from U(R X ) = 0 we arrive at 
R x ss ( 3 _ 1 C 3 n 8 ) 1 / 3 . For n = 100 this estimate with 
C 3 = 1 leads to R x ss 8 pm in a reasonable agreement 
with our computed value. Further, we evaluated molecu- 



g + r 


FIG. 2: A laser pulse nominally resonant with the r — g tran¬ 
sition can resonantly excite molecular states. The molecular 
wave-packet is efficiently excited within a window A R x deter¬ 
mined by the slope of molecular potential U(R) and molecu¬ 
lar Rabi frequency f2 m . Once excited the wave-packet rapidly 
accelerates out of resonance and rolls down the slope of molec¬ 
ular potential. 

lar Rabi frequencies = £, m ^o (typically a fraction 
of Ho, see Table [i]). from the eigenvectors of the numeri¬ 
cally diagonalized molecular Hamiltonian. 

The atoms are efficiently laser-coupled to the molec¬ 
ular resonances only in a small window of R , when the 
detuning U ( R ) is comparable to H m (see Fig.[2j). Thereby 
we define an effective radial width of the molecular reso¬ 
nance 

AR X =n m (R x )/\U\R x )\, ( 1 ) 

where S2 m and the derivative of the molecular poten¬ 
tial U'(R X ) = dU(R x )/dR are evaluated at the reso¬ 
nance crossing. For the outermost resonance A R x ss 
HoUn 20 / 3 (c 3 /3) 1 / 3 /9. Each molecular resonance there¬ 
fore defines a “resonance shell”, a spherical shell of radius 
R x and radial width A R x centered at a given atom. The 
average number of atoms inside the resonance shell, 

A N x = 4TrR x AR x pd , (2) 

is a relatively small number in a typical experiment. For 
parameters of Table [TJ the outermost resonance shell con¬ 
tains less than ~ 0.1 atoms. 
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The presence of molecular resonances implies several 
consequences for the Rydberg blockade, the most im¬ 
portant being the atomic loss. Indeed, inside the shell 
two excited Rydberg atoms are a subject to a mechan¬ 
ical force —U'(R X ). This force can be either attractive 
or repulsive. (An example of the repulsive resonance is 
at R x ~ 6.2 pm, see Fig. [2]). The diatom would sep¬ 
arate into two 99s and I01s Ry atoms with a kinetic 
energy of relative motion equal to the dissociation limit 
S m ~ — 3?r -4 , which is ~ 10mK for n = 100. These 
atoms may escape the trapping volume, effectively re¬ 
ducing the number of blockaded atoms. Since the atoms 
are accelerated out of the resonance shell on timescales 
T a = ^AR^M a /\U\R x )\ « 1/fio, we adopt a sim¬ 
ple model that once a pair of atoms is promoted to a 
molecule, the associated motional wave packet quickly 
leaves the resonance shell with its atomic constituents no 
longer interacting with the laser field. 

Now we analyze the dynamics of collective atomic en¬ 
semble excitations. We will derive the expressions for 
the damping (loss) rate in two approximations: (i) as¬ 
suming that atoms are frozen in space (static limit) and 
(ii) collisional model (impact limit). Remarkably both 
approximations yield identical loss rate. We consider an 
ensemble of N atoms initially in the collective ground 
state |G) = |gi <72 ■ • • <7jv)- A laser pulse couples |G) to 
a superposition of singly excited Rydberg atoms \R t ) = 
|< 7 i • • • gi_irigi+i ■ ■ ■ gw}- These atoms can be further pro¬ 
moted to doubly-Ry-excited diatom states M l3 ), involv¬ 
ing atoms i and j. There are N m = N(N — l)/2 diatom 
states, with their resonance detunings A (Rij) = U(Rij) 
and Rabi frequencies ( Rij ) determined by their inter¬ 
atomic separations Rij. Expanding the total wavefunc- 
tion in this basis (wo is the laser frequency resonant with 
the g — r transition) 

N N N 

I*) = c g e iu>ot \G) + ^ a\Ri) + e~ iu>ot £ £ m, 3 |Afy> 

i i j>i 

and applying the rotating-wave approximation, we arrive 
at 
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When all molecular detunings A^ are large, the sys¬ 
tem undergoes the ideal Rabi flopping between the col¬ 
lective ground state and symmetric combination of sin¬ 
gle Rydberg excitations ( Cg(t) = cos(v / AAIoi/2) and 
c\(t) = —ism(\fNSlnt/2)/yfN). We focus on the av¬ 
eraged collective dynamics and introduce collective am¬ 


plitude c s so that d(t) = c s {t)/\/N. We assume that 
all ffd = owing to the weak R. l3 dependence inside 
the resonance shell. The collective amplitudes satisfy 

(fiff = VN£l 0 ) 


ic a = -Co . 

9 2 

.. 

*c s — 2 C 9 


- - 



inn. 
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2VN 
)rriij + (f l m /y/N)c, 


( 3 ) 


Now we fix the positions of all atoms (this require¬ 
ment is relaxed later on) and split the time axis into 
time intervals consisting of a short laser pulse of du¬ 
ration t p < 1/ft* and field-free acceleration time r a 
during which the excited diatom wave packet leaves the 
shell. Because of the mechanical forces the molecular 
amplitudes inside the shell are reset to zero values be¬ 
fore the next pulse arrives (this is reminiscent of the 
Markov approximation [TS]1. By taking the limit r a —>■ 0 
we arrive at a continuous Rabi drive. Integrating the 
last equation over time interval (t, t + r p ), one obtains, 
rriij(t + T p ) = c s (t)(fl m /VN) {exp(iAijT p ) - 1} /Aij. 
We have set m^(t) = 0 as discussed. Notice that 
the r.h.s. spikes at A^ = 0, i.e., within the reso¬ 
nance shell. Ensemble averaging yields ( rriij(t + St)) = 
-ic s (t)n(Q m /y/N)(4:irR 2 x )/(\U'(R X )\V S ), where V s is the 
blockaded ensemble volume. By substituting this rela¬ 
tion into the equation for c s we arrive at a set of damped 
equations (non-Hermitian Schrodinger equation) 


n N 

Us = C g - i'tm.Cs (4) 

with the molecular-resonance (amplitude) loss rate 

7m = 7rAA x n m /2 = 2n 2 pdCnRl^l/W{R x )\ . (5) 

For the outermost resonance, 7m ~ 2t 2 pd^n^i^ n 12 /27. 
The above derivation neglected atomic motion and is 
valid for very cold ensembles. In the SM, we take into 
account the thermal motion of the atoms, using the im¬ 
pact approximation. We find the result for 7m that is 
identical to Eq. 

In addition to the loss, the same f1 drive induces an 
AC Stark shift which is different for states in which dif¬ 
ferent atoms are excited. This inhomogeneous broaden¬ 
ing results in an additional loss of coherence of the Rabi 
oscillation. In the SM, we show that this broadening is 
given by 

5A^(7n/2)p d e m R 2 x n 2 /\U\R x )\ (6) 


in the limit where R x is much smaller than the size of 
the ensemble. Although the dependence of the system 
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parameters is the same as for 7 m in Eq. ([5| , the prefactor 
makes 6 A roughly two times smaller than j m . During the 
Rabi dynamics, the broadening makes the signal decay 
as exp[— (fJA) 2 ] over time. 

The rate formula is to be summed over all resonance 
shells: 7 ^ =_V fc 7 m) where 7 ^ is the individual shell 
contribution ([ 5 JT For our example in Table [ij 7 ^* is en¬ 
tirely dominated by the outermost crossing. The rea¬ 
son for this prominence is that at smaller R. the poten¬ 
tials become steeper and the molecular Rabi frequencies 
become diluted. Notice that the long-range molecular 
Hamiltonian used in computing the molecular potential 
curves in Fig. [I] holds only for R 2n 2 ao ~ 1 pm for 
n = 100 , i.e. when the electronic densities do not overlap. 
However, all the qualitative arguments that the molecu¬ 
lar excitation rates should be suppressed compared to the 
outermost resonance shell still hold even for small R. 

The rate scales steeply with n, "f m oc n 12 . In particu¬ 
lar, it is commonly believed that the blockade fidelity can 
be improved by going to high-n Ry states, because the 
probability of off-resonant Ry excitations is suppressed 
as n ” 22 in the van der Waals blockade. We see that in¬ 
creasing n while suppressing off-resonant Ry excitations 
also increases the undesired molecular loss rates. 

Eqs. © reduce to the damped oscillator equation of 
motion c s + (fl// /2) 2 c s — 7 m c s = 0 with solutions 

c s (t) = -i (f iff/fl d ) sin (fl d t/2) , (7) 

c g (t) = (cos ( fldt/2) + ( 7 m/fld) sin {fl d t/2)) e ~ J ^ t , ( 8 ) 

where fl d = Hq (l — rj 2 ) , 77 = 7 m /fl$ ■ Driven ensem¬ 
ble exhibits damped collective Rabi oscillations with a 
frequency fl d < Q (J V . One may distinguish between three 
classes of solutions US: under-damped (77 < 1 ), critical 
(77 = 0) and over-damped (77 > 1). Explicitly, 

77 = 2 * \ s mR x )\ ^ R * n °« nl2 ^°- 

Thereby increasing Do, n, or p d can cause the ensemble 
to exhibit over-damping of collective Rabi oscillations, 
at which point they no longer resemble oscillations. In 
the under-damped regime, the loss per collective Rabi 
cycle determines collective qubit operation fidelity F = 

(flo /fl d ) 2 e~ 2mi ■ 

The molecular loss can account for some experimen¬ 
tally observed imperfections. E.g., Dudin et al. [TO] 
have effectively measured the damping constant for col¬ 
lective Rabi oscillations in a mesoscopic ensemble of 102s 
Rb atoms. They found that the Rabi oscillation loses 
10 — 20% of its contrast in a single cycle. Our calculation 
can account for a loss of ~ 5% during a single oscillation 
cycle. While the agreement seems to be adequate, we 
emphasize that it may be fortuitous as the experiment 
has been carried out in the presence of magnetic field 
which was excluded in our analysis and would introduce 


additional resonances. We also neglected the s — d ex¬ 
citation channels (allowed in the excitation scheme [ 10 ]) 
when computing molecular Rabi frequencies. Moreover, 
the experiment m is affected by a multitude of other 
decoherence effects. At this point it may be desirable to 
design experiments that could disentangle various deco¬ 
herence mechanisms, and the molecular losses in partic¬ 
ular. 

While the total number of atoms N remains constant 
during the coherent evolution, the wave function acquires 
out-coupled diatom wave packets. If the measurement of 
the total number of atoms were to be made, the number 
of atoms remaining in the ensemble would be N e ff(t) = 
N{\c g {t)\ 2 + |c s (f)| 2 ). By manipulating Eq. 0, one 
finds that N e g = — 27 m |c s (t)| 2 fV, resulting in N e g(t) = 

* {i + 2 (?£) 2 sin 2 (^) + (^) sin(D d t) J or 

averaging over many cycles N e s{t) = IV (1 + 
ilm/fld) )e -7mt , i.e., the effective number of atoms re¬ 
maining in the ensemble decays exponentially. The 
quantity 1 - (|c g | 2 + |c s | 2 ) = (N - N eS (t))/N ~ 
('ym/fld) 2 e~ lmt also determines “leakage” from the col¬ 
lective qubit space. Clearly, to minimize the leakage one 
has to require that 7 m t <C 1 or 7 m <C floy/N. For pa¬ 
rameters of Table [i] the coherent evolution is limited to 
t < 10//S. 

One may visualize the “leakage” from the collective 
qubit space as a modulated outflow of molecular wave 
packets from the blockaded volume. As an illustration, 
the outermost molecular resonance produces admixtures 
of 101s and 99s Rydberg atoms. If the ensemble is 
trapped, the out-coupled (di)atoms may linger inside the 
ensemble depending on the released kinetic energy and 
the trapping potential height. Such atoms do not res¬ 
onantly interact with the laser field of the Rabi drive. 
However, they do interact with the remaining ensemble 
leading to energy shifts through the interactions with the 
remaining atoms. Such mechanisms can be also relevant 
for untrapped ensembles, where the outflowing diatom 
wave packets may interact with the remaining ensemble 
while transiting out through the volume. In addition, the 
present discussion focused on Rydberg S-states, the un¬ 
desired effects can be enhanced for Rydberg states with 
higher angular momentum. This is due to the presence of 
closely spaced states, due to, e.g., spin-orbit interaction, 
that can result in molecular crossing at larger R. 

The unwanted effect of molecular resonances can be 
suppressed by using tight traps for individual atoms prior 
to excitation, such as optical lattices. The idea is po¬ 
sitioning atoms such that excitation to molecular reso¬ 
nances is not allowed. As shown in the SM, the loss can 
be suppressed by a factor of ~ 100 , if the tightly trapped 
(~ 20 nm) individual atoms are prepared in a 3D optical 
lattice with the lattice constant tuned to avoid the reso¬ 
nances. By choosing the lattice constant, the outermost 
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resonant shell R x = 6.2 /jm can fall in a gap between 
density peaks, largely reducing its effect. High fidelity 
manipulation of Rydberg atoms in a lattice has been ob¬ 
served recently ns], 

We investigated how molecular resonances limit the fi¬ 
delity of Rydberg excitations in an atomic cloud. Under 
continuous driving pairs of atoms can be promoted into a 
doubly excited Rydberg states, if they are separated by 
certain resonant distances. These resonant pairs repel 
each other and leave the cloud. To mitigate this detri¬ 
mental effect, trapping the atoms in a tight optical lattice 
can be used, where they are kept away from resonance. 
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Supplementary Materials 


DERIVATION OF THE MOLECULAR EXCITATION RATE IN THE IMPACT APPROXIMATION 

As discussed in the main text, the collisions leading to strong coupling to molecular states are short and well- 
separated. Let’s consider one of such collisions of atom i with an atom j. The molecular probability amplitude 
satisfies the equation 


irhij = Aijiriij + (q + cj). 


( 9 ) 


The detunings A^ (molecular potentials with the zero energy at the |r) + |r) dissociation limit, A *j = U (Rij (t))) 
and Rabi frequencies are time dependent because of the atomic motion. Introducing 


dt' ) , 


we recast Eq.([9]) into 


rriij ( t ) = rhij ( t ) exp ^-i J (t 7 ) 


*-tt mij (t) = (cz + cj) exp [i A tj (t 1 ) dt 


“dt 10 w 2" 

We are interested in the probability of molecular excitation due to a single collision, 


( 10 ) 


P m = \rriij (oo)| 2 = |rhij (oo)| 2 . 


Integrating Eq.(10), we arrive at 


mi 


(oo) = —i 


2 n m (t) (D {t) + Cj (t)) exp ( i 


U (Rij (t')) dt' ) dt 


( 11 ) 


To evaluate this probability we approximate Rij ( t ) with straight-line trajectories, 

= {(u {t - t c )) 2 + p 2 } 1/2 . 

Here p is the conventional impact parameter, v is the relative atomic velocity, and t c is the time of the closest approach. 
The atoms reach the resonance region when Rij(t) = R x . Clearly, one has to require that p < R x for this to occur. 
The associated moment of time t x is 


c ) = ±\JR 2 X -p 2 - 
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The exponential in the integral e rapidly oscillates except when the phase U {R l: j it')) dt' is stationary. The 
prefactor varies slowly in time compared to the exponent. This forms the basis for evaluating using the stationary- 
phase method. Let’s review the basics of this method. Consider an integral 

/ +oo 

9(t) e^dt, 

-OO 

where g (t) varies slowly compared to the rapidly oscillating exponent. The main value of the integral is accumulate 
in the regions where the phase is stationary, i.e. 


dcj) (f*) /dt = 0. 


Expanding the phase in the vicinity of the t* 

</> (f) ~ 4> (4) + ~4>" (4) (■t - 4 ) 2 . 

Then 

I « g (4) J exp (4) (t - 4) 2 ^ dt = 

9 ^ exp OJ si9n ^ (**))) 

In our case cj) (t) = U ( R t j (t')) dt', and the stationary points correspond to the crossing of the resonance shell 

dcj> (4) /dt = U ( Rij (4)) = 0, 

i.e. 4 = t x . Notice that we have two stationary points corresponding to two crossings of the resonance shell. The 
two times are separated by 


tl - f x = 


‘hlK - p 2 


In general, both points can contribute. However, once promoted to the molecular state, the atoms experience strong 
mechanical forces and are accelerated out of the resonance. Therefore we will neglect the interference effects when 
computing the probability p m and add the two contributions incoherently (this provides the upper limit on P m ) 


P m = 7 \a (t x ) + Cj {t x )\“ 


1 


W ( 4 ) 1 , 

Further we evaluate the second derivative of the phase evaluated at crossing points 

dU dRjj 


d (j) (t x ) /dt 2 = 

Urn V 


dRtj dt 


A R x Rs, 


1 (' t x - t c ) = 


A R x R x 

From here one could define the effective duration of collision 


R\ - p 2 


A R x 1 




or 


Pm — 7T \Ci ( t x ) + Cj (^ai)| 


A Ry 




v 


Rl - p 2 
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We further approximate the time evolution of single-Ry-excitations via their un-coupled time evolution, c(,(f) = 
i/y/N sm(y/Nfl 0 t/2)) 


47 r 

Pm ( P > t ) = — fir 


A Ry 


Rx 


\Jr\- p‘ 


: sin 2 (VNn 0 t/2)). 


p m = 0 for p > R x in the stationary-phase approximation as there are no crossing through the resonance region for 
such impact parameters. 

Number of atoms lost due to a single collision is AN = —2P m . Now we sum the probabilities over multiple 
collisions. Number of atoms in a relative velocity group dv passing through the area 2tt pdp per time interval dt 
is equal to 2nnpdp\v\ f (v) d 3 vdt, where n is the number density and f (v) is the velocity distribution. Then the 
compound atom loss satisfies the equation (here the factor of 1/2 is introduced to correct for double-counting) 


dN „ . , N . , Ar 

— = -2 7m (t) — = -7m (t) N, 

7 m(t) = J J 2ttup dp \v\ f (v) d 3 v p m 


Explicit evaluation yields the cross-section 

&m (t) = 2n f p dp p m (p,t) = (27t R 2 x ) sin 2 (v / ]vnot/2)). 

Jo N v 

and the rate 

7 m ( t) = n (27r R 2 X ) —fl m AR x sin 2 (VNn 0 t /2)) 

For a spherical volume of radius R s , n = 3 N/ ( 47 ri?f), thereby 


7 m (f) = <S» (^) 


sin 2 (y/NSlot/ 2)) 


The rate equation has the solution 


N(t) = N ( 0 ) exp J 7 m {t') dt'^j 

in NQot 


— Tm 1 

r 

= 37r 

(Rx 
V Rs 


2 y/Nflo 




For sufficiently long time (t 47 t/ ^y/NQ,oj), the total number of atoms falls off exponentially as 

N (t) = N (0) exp (- 7 m t) ■ 


Finally, the experiments are carried out with mesoscopic ensembles and as discussed in the main text, the radius of 
the resonance shell R x maybe comparable to R s (or blockade radius). It is clear that if R x > 2R s , the atoms are not 
going to be affected by that particular molecular resonance. In analogy with Eq.() of the main text we may further 
introduce a geometric probability factor g ( R x /R s ). 

7m (t) 7m (t) g (Rx/R s ) 


Further rates from multiple resonance add 

7m (t) W 

k 


where 7 ^ is the rate due an individual resonance shell at R x . 



DERIVATION OF THE INHOMOGENEOUS BROADENING 


Hamiltonian 

We assume that the N identical atoms move negligibly over the entire extent of the dynamics in question, so we 
need to track only the electronic degrees of freedom. We model each atom as a four-level system, with states |g), 


|r), | r'), | r"). Let us define the following collective states, 

|G) = | g)® N , (12) 

I j) = crj|G), where <7j = \ 9 )j{r\j, j = 1,2,... iV, (13) 

\j,k) = where ct' = \g)j{r'\j, and a” = \g)j(r"\j, j^k (14) 

An external driving field coherently couples | g) with |r), |r') and | r"). When two atoms are in | r'r") or | r"r f ) states, 

they interact via the Rydberg interaction. The resulting Hamiltonian is 

\ H = Y -^(l G )0'l + h - c -) + X] yp(U>0'j fe l + \j)(k,j\+ h.c.) +J2 A jk\j,k)(j,k\ (15) 

j j,k j,k 

In the ideal case, when A —> oo, |G) is coherently coupled to the symmetric combination of a single excitation, 

= (“) 


and the resulting dynamics is a Rabi oscillation between |G) and IS 1 ), with Rabi frequency Qr = y/NQ, if the system 
starts in |G). 

To investigate the deviation of the real dynamics from the ideal one, we focus on the coupling of | j, k) states to \S). 
Using this notation the Hamiltonian can be written as 

\ H = ^^(|G}<S|+h.c.) + Y, %(|5)<M ifc |+h.c.)+ £ A ok\Mj k ){M jk \, (17) 

j,k>j V j,k>j 

where | Mjk) = \k k )+\k,j) ^ an£ j the non-symmetric combinations are not coupled to jS 1 ). 


Broadening 


We adiabatically eliminate the doubly excited states {\j,k)} from Eq. to arrive to the effective Hamiltonian, 


n 


j,k>j 


^H c ff = hn 0 + Y[ A j - ) \ j)(j\ + Y ^(U)(*l + l fc )(j'l) 


(18) 


where hn 0 is the first term in Eq. (|17[), and the new coefficients are 


A, 

'P 

h 4 

(19) 


n 2 


£ljk 

2A j k ’ 

(20) 

r, 

= 2tt y ( Q m 6(Ajk), 

(21) 


k^j 


where V indicates principal value. The Tj terms describe the resonant excitation, which is as discussed in the previous 
section. 

The j dependence of A j results in inhomogeneous broadening, <5A := y/ (A 2 ) — (A) 2 . 


o 2 

SA = -^ 

4 


\ 


Var, £ _ ^ l = 


M4 

4 


\ 


Var ( Y 
k^j 


V 


R % 


A E ft'.': - St 


jk 


( 22 ) 
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where we used hAjk = U(Rjk) = jj£ — A e, he. a van-der-Waals interaction potential between the Rydberg atoms, 

jfc 

and eliminated Cg by using that R x = . 

The sum can be written as 


D:> ■= E 


V 




a e bSu - Ri 


jk 


N [« N f ,3 Vr 6 

VA E J rfc |r fc - r^l 6 - R 6 X VA E J r r 6 -i?®’ 

v v ' 


(23) 


where the integrand is more conveniently written in terms of r = |r| = |r^ — r ; |. r can be seen as ‘local’ spherical 
coordinate, centered around rj. Due to global rotational invariance of the problem we can set r j = rj z without loss 
of generality. We then find (see appendix [| 


DjV A e 
N 



o 


Vr 6 

r 6 - R 6 X 


R+r d 



R—Vn 


Vr 6 

r e - R 6 x 


r 2 + r? — R 2 \ 
-- dr. 

2 rrj j 


(24) 


We now assume that the singularity is in the first integral, R x < R — rj. As a result, the second integral is no longer a 
principal value integral and since r > R — rj > R x we will furthermore approximate « 1. The second integral 

is now straightforward to evaluate. Finally we use the indefinite integral 


1 


1 


r 6 — Rt 


dr = -r 3 — -R^arctanh ( 


to find the remaining principal value integral. The result is 


N 


A e 


1 - 



3 

arctanh 




(25) 


(26) 


Now, we can determine the averages (D) and ( D 2 ), but since the above expression for Dj is only valid when rj < R—R x 
we modify the averages to only average over a sphere with radius R — R x : 


(D) 
( D 2 ) 


1 

N 





V 


3 

2 (R-R x ) 3 


R—Rx 



0 



3 

2 (R-R x ) 3 


R—R x 

J r j{Dj) 2 drj, 

o 


where we also used that Dj only depends on Vj = r ;/ 1. Using the expansion 


D j = . 

J A 


N 

E 


1 - 


R% 


R 3 (R — rj ) 3 


+ ... 


up to the first term including r, the averages are found to be 


, , N 

- Ae 

m = (A 


3 (R* 
2 l R 


_l_ 3 ( R x 


R 


R 


R 


1- 3 +3 


It, 


R 


where we made an expansion in powers of As a result we find (lowest order in ^) 


SD = V(-D 2 > - ( D ) 2 


N fz ( R x \ 7/2 

A e V 5 V R ) 


(27) 


(28) 


(29) 


(30) 

(31) 


(32) 
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With this result we can write the broadening 5A as 

<33) 

where the power of 5/2 was also reproduced using numerical calculations. If we apply the same calculations on a 
dipolar potential, U(Rjk) = 4 ? — A E , the result is 

;?'fc 


SA 


-nhClz 


,Pd 


Ri 

I U'(R X )\' 


which surprisingly has the same dependence on R x as 7 m . 


(34) 


Effect on coherence 

The effect of this inhomogeneous broadening is well approximated by an additional decay of the coherence between 
the ground state |G) and the symmetric single-excitation state | S) = I Ri)/N, by a factor of e~( SAt ) / 2 . This is 
a much weaker effect than the pure exponential decay, set by j m , and since SA < -y m for the parameter regime in 
consideration, the effect of inhomogeneous broadening can be neglected as long as 7 m t < 1 . 


DEPHASING 


In the absence of molecular coupling, Rabi flopping occurs between the collective ground state and the collective 
single Rydberg excitation 

i# 7 (*)> = Y c *( t ) I-* 2 *) = isin f Y2 i^) 

k—1 \ / V k—1 

\& (t))=cl(t)\G)= cos \G) 

\* I (t)) = \G I (t)) + \R I (t)) 

With the molecular resonances 

N 

i-R o)) = Y2 ( c fe w+ ak (*)) 1 ^) 

k=1 

I G(t)) = (cl C t) + a g (t)) \G) 

\*(t)) = \G(t)) + \R{t)) 

where we introduced corrections ak (t) to the ideal amplitudes cj, ( t ). In the ideal case 
The overlap between the target | R 1 [t)) and the actual | R (t)) states reads 

N 

(T J (f) IT (f)> = (cl ( t))* (cl (t) + a g (t)) + Y (4 {t))* (cl ( t ) + a k (t)) = 

k=l 

N N 

= ( c 9 (t))* ( C fl 0) + a 9 0)) + Y l C fe (0| 2 + Y ( C fc W)* ak (*) = 

fe=l 

N 

= 1 + (4 ( 0 )* «9 (*) + (*))* afc (0 

fc=l 
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Distance between |'J> / (t)) and |\E r (f)) 

d 2 = (V - 'I' 7 |’L - ^L 7 ) = 1 + (T (t) 1$ (t)) - 2Re(<L 7 (t) (t)) = 

l N 

((^O) 1^ (*)) - 1) - 2 Re (c^ (f))*a s (t) + ^ (4 (t))*a fc (t) 

\ k=1 


(0 1^ (*)> = 1 + 2Re | (c 7 (t))*o fl (i) + ^ (4 W)* a fc 0)) + |a s | 2 + ^ |a fc | 2 

\ fc= i / fc 

^ 2 = m 2 +X! m 2 

k 

As in the derivation of the molecular exctition probability, consider a transient coupling to a molecular state M t j). 
Consider time evolution of the rydberg amplitude of atom z. 

iC-i —QftCg T —Qrn'Tflij 
1 N 

idg 2 ^^ ^ ’ 

irhij = Aijniij + ( Ci + cj). 

or with Ci = c\ + a*, we find (here we take into account that the time evolution of and a j are identical - this is true 

for the first collision only) 


as derived in Sec. 


XCLg — 

ibji = — + —Q-rnTYlij 

XTflij — + Q'ffiCj 


ij ( t ) « rny (f_) — ittmcj ( t ) exp ( t )) J exp (z</> (t')) dt' 


rriij (t-) - iflmcj (t) exp (-i<£ (t)) 


Since 


d 2 Id Id 
^ — 2^° dt 2 m - - ^ 


dt 2 

d 2 


dt 


*7 


dt 2 


I want to neglect Ay-my. Let’s do this first 


. 1 rV 2 1_ d 

z— TrOi = —z-i2 n aj H—s2 m — mi,- 

dt 2 2 0 2 dt 3 

CLi -j- —flQnj = Z —(A^jZTZjg + D m C,^ 


d 2 1~ 2 1 

dt 2 " n+ 2 r ° a,: - 2 U 


/ v 7 NQ(jt \ 


we can consider the two contributions separately as the solutions linearize in the inhomogeneous part anyways. 
Basically these oscillations persist even without the collisions (I guess one needs to look at time-dependence of f2 m ). 
The second part 

d 2 l^o l n * 

^^ 2 ^ 1 A 2^o a * — 2 znjy 
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Actually let’s leave the eqn in this form (to use integration by parts) 


d 2 1 „ 2 .1 „ d 

d^ a ' + -n 0 Oi - -*2 n ">dis¬ 


solution wo = fio/v 7 ^ 


(t x + St x ) = b+e iuiot + b-e~ iuot + 


i /n 


4 V w o 

i /n 


o -iuj 0 t+ 


r t x 


exp (+iuJot') to (t') dt' 


4 V w 0 


™ | e +i«ot+ 


t x Si x 
rt x -\-St x 

f t x —St x 


exp (—iiijot 1 ) rh ( t') dt' 


The initial conditions are determined by 


x i ( t x - 5t x ) = b + e lUot - + fo_e *“ ot - 


f 

f*t x -\-St x 


rt x +<5^x 


exp (+«wo<) to (t 7 ) dt 7 = exp (+iuiot) m (i)| t + — *wo 


exp (+iwot 7 ) to (t 7 ) dt 1 


t x St x 


-ZUJo 




exp (+iu>ot) m (t)|*+ = 
exp (+iwot 7 ) to (t 7 ) dt' = —iojorriij (t_) 


- (-iw 0 ) iQmc{ ( t ) 


/»ta;+<5t x 

I exp (+iwot 7 ) dt'+ 

‘t x ~\~&tx 


r tj St x 


exp (+iojot) exp (— i<j> (t)) dt / exp {i<j> (t 7 )) dt' 


EFFECT OF ATTRACTIVE POTENTIAL 

The attractive potentials (see at R x ~ 4.6 /im in Fig. 1 of the paper) can lead to auto-ionization in the small R 
region nzi- Such a process would free an electron and a molecular ion, with their Coulomb fields blockading the entire 
sample. While in our illustrative example, laser coupling to attractive potentials is negligible, it may be not be the 
case in general. Qualitatively, to reduce auto-ionization one needs to pick Rydberg states such that the potentials 
inside the most strongly coupled resonance shells are repulsive. 
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SUPPRESSION OF RESONANCE BY TRAPPING IN REGULAR LATTICE 


We investigate the possibility of trapping the atoms in an 3D cubic optical lattice, in order to suppress the effective 
decay rate 7 , due to double Rydberg resonances. We use the homogeneous atom density case as a benchmark, 

Thom = J^Tf = 27r 2 n 0 p^^i?x,iAi?x,j, (35) 

3 » 

where the summation is performed over all resonances at the crossing distances R x ,i , each having a width of A R x ,i 
and molecular coupling factor Here p is assumed to be constant. 

We assume that in a 3D cubic lattice, each lattice site holds a single atom, trapped in the ground state of the 
harmonic trap. Let a denote the lattice constant and d be the size of each trapped wavefunction. The deeper we make 
the lattice, the smaller the d/a ratio can be. For given d, a values, we can plot the 3D density p(R) of the atoms as a 
function of distance from a particular lattice site. This is shown on Fig. [3] with a blue curve. 
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0.44 

0.13 


0.091 0.55 



FIG. 3: 3D density of the atoms in a cubic lattice (lattice constant: a = 0.995 pm) as a function of distance, R. Each atom is 
confined by a harmonic trap to a region of size d = 0.02 pm. Red vertical lines indicate the position of the resonances given in 
Table 1 of the paper. The numbers shown next under the lines are the moleculare coupling coefficients, for each resonance. 

By taking the R dependence of the atom density p(R) into account, we can write the total decay rate as 

Tlattice = ^Tf = 27T 2 f2 o y^p(J? X| i)^i? 2 ,j&Rx,i, (36) 

3 i 

numerically evaluate, and compare it with the homogeneous result, 7hom- On Fig. 0 we plot Tlattice/Thom as a function 
of the trap confinement d for different fixed values of a. A confinement of d/a ~ 0.01 can suppress the decay to 2-10% 
of its homogeneous density value, depending on the accuracy of the fine tuning of the lattice constant. 
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cd 


<F~ 



d (/im) 

FIG. 4: Decay rate in a 3D lattice normalized with the homogeneous decay rate result for lattice constants a = 1.00, 0.995 
and 0.99 /im. The suppression is strong for small trap size, and diminishes (7i a ttice —> 7hom) for large size. Small changes in a 
result in significant changes in the suppression. This is due to the detailed peak structure of the 3D density in the lattice. As 
a result, fine tuning of the lattice constant is required. 


Local spherical coordinates 


Since we make the substitution r = r*, — rj the boundaries of the polar integral change. Depending on the relation 
between r, r, and R there are three regions, as shown on Fig. [5j First, if r < R — 77 , then the entire sphere of radius 
r = |r| lies within the boundaries of the cloud, and therefore we have 0 < 9 < ir. Similarly, if r > R + Tj , then the 
opposite is true: the cloud lies entirely inside the sphere of radius r, and therefore there is no contribution to the 
integral from this part. 



r < R - rj 


R — rj < r < R + rj 


R + r-j < r 


FIG. 5: We only need to integrate over that part (red) of the dotted sphere that is inside the cloud. This corresponds to a 
modification of the boundaries for the 0-integral, depending on the relation between r, rj and R. 


If R — rj < r < R + r-j, then there exist a circle, where the sphere of radius r intersects the boundary of the cloud, 
as shown on Fig. [5} The angle between the segments rj and r is 


/ r 2 + - R 2 \ 

^ 2 ^ ) ’ 


7 t ~ 6 C = arccos 


(37) 
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and therefore we have 9 C < 9 < n. Thus we can rewrite the integral in Eq. (231 as 


R—rj 


2-rr 


R+rj 


2 t r 


Jd 3 4 5 6 7 8 r f(r) = J r 2 dr J sin (9)d6 J d(j)f(r)+ j r 2 dr J sin (9)d9 J d(/>f(r) (38) 

V' 0 0 0 R-rj e c 0 

where f(r) is the integrand that only depends on r. Evaluating the integrals over <f> and 9 we find 

R—r-j R-\-rj 


lb i j “T'j / 

Jd 3 r f(r) = J 47 rr 2 f(r)dr+ J 2nr 2 f(r) ( 1 — 

V' 0 R-r; ' 


r 2 + r" 2 — R 2 \ 

- w 1 - \dr. 

2 rrj I 


(39) 
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